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Abstract 

A general way to construct ladder models with certain Lie algebraic 
or quantum Lie algebraic symmetries is presented. These symmetric 
models give rise to series of integrable systems. It is shown that cor- 
responding to these SU{2) symmetric integrable ladder models there 
are exactly solvable stationary discrete-time (resp. continuous-time) 
Markov processes with transition matrices (resp. intensity matrices) 
having spectra which coincide with the ones of the corresponding inte- 
grable models. 
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Integrable models have played significant roles in statistical and condensed matter 
physics. Some of them have been obtained and investigated using an algebraic or coor- 
dinate "Bethe Ansatz method" |jl], The intrinsic symmetry of these integrable chain 
models plays an essential role in finding complete sets of eigenstates of the systems. On the 
other hand, stochastic models like stochastic reaction-diffusion models, models describ- 
ing coagulation/decoagulation, birth/death processes, pair-creation/pair-annihilation of 
molecules on a chain, have attracted considerable interest due to their importance in many 
physical, chemical and biological processes [0. The theoretical description of stochastic 
reaction-diffusion systems is given by the "master equation" which describes the time 
evolution of the probability distribution function This equation has the form of 

a heat equation with potential (i.e., a Schrodinger equation with "imaginary time"). If 
an integrable system with open boundary condition can be transformed into a stochas- 
tic react ion- diffusion system, e.g., by a unitary transformation between their respective 
Hamiltonians, looked upon as self-adjoint operators acting in the respective Hilbert spaces, 
then the stochastic model so obtained is exactly solvable with the same energy spectrum 
as the one of the integrable system @, 

In [Q, we have presented a general procedure to construct open chain models having a 
certain Lie algebra or quantum Lie algebra symmetry by using the coproduct properties of 
bi-algebras. These models can be reduced to integrable ones via a detailed representation 
of the symmetry algebras involved. In recent years spin ladders have attracted consid- 
erable attention, due to the developing experimental results on ladder materials and the 
hope to get some insight into the physics of metal-oxide superconductors 0. In this 
letter we study the construction of ladder models with certain Lie algebra or quantum 
Lie algebra symmetry. We show that the integrable quantum spin ladder model dis- 
cussed in |T0| can be obtained in this way and it can be transformed into both stationary 
discrete-time (discrete reaction-diffusion models) and stationary continuous-time Markov 
processes with transition matrices resp. intensity matrices having the same spectra as the 
ones of this SU{2) invariant integrable ladder model. 

Let A be an associative Lie bi-algebra with basis e = {e^}, a = 1,2, ...,n, satisfying 
the Lie commutation relations [e^, C/j] = C^/jC-^, where C^p are the structure constants 
with respect to the base e. Let A (resp. C(e)) be the coproduct operator (resp. Casimir 
operator) of the algebra A. The coproduct operator action on the Lie algebra elements is 
given by Ae^ = Cq, (S> 1 + 1 (8> ea, 1 stands for the identity operator. It can be immediately 
checked that [Acq, Ae/3] = C^^Ae-y and [AC(e), Acq,] = 0, a = 1,2, ...,n. 

Let us consider a two-leg ladder with L rungs. To each point at the i-th rung, i = 



2 



1, and ^-th leg, 6 = 1,2, of the ladder we associate a (finite dimensional complex) 
Hilbert space Hf. We can then associate to the whole ladder the tensor product Hi 
Hf C?> H2 (S> i/| ... (S> Hj^ Hi. The generators of the algebra A acting on this Hilbert 
space associated with the above ladder are given by = A^^'^Cq,, a = 1,2, ...,n, where 
we have defined 

A"" = ( 1 . . . O 1 ^A)...(l 01^ A)(l 8) A)A, MmelN. (1) 

m times 

Ea also generates the Lie algebra A: [E^, E/^] = C^pE^. 
Let 

3 2 

/i = 5:^a,A'A]AC(e), (2) 

i=i j=i 

where A^ = (A®1), A^ = (lOA), A^ = (AolOl), A^ = (1® A®1), A^ = (lOloA), 
G (C such that h is hermitian. Let IF denote a real entire function defined on the 2L-th 
tensor space A A (^^ ... A of the algebra A. We call 

L-l 

H=Y^ lF{h),,,+, (3) 

i=l 

the (quantum mechanics) Hamiltonian associated with the ladder. Here lF{h)i i_^.i means 
that the four-fold tensor element lF{h) is associated with the i and i + 1-th rungs of the 
ladder and acts on the space H^ ® Hf Hl_^-^ (g) Hf^^, i.e., 

lF{h),^,+i = Ij (g) 1? ® ... ® o o iF(/i) o 11^2 ® l'+2 ® ••• ® li ® li- (4) 

[Theorem 1]. The Hamiltonian if is a self-adjoint operator acting in Hi ® Hf ® H2 ® 

i/| (g ... (g) Hi (g) and is invariant under the algebra A. 

[Proof]. That H is self-adjoint is immediate from the definition. To prove the invari- 
ance of H it suffices to prove [H, E^] = 0, a = 1,2, ...,n. 

From the formula for the above coproduct we have 

L 
i=l 

where (e„)i = 1^ ® 1^ o ... o o if^^ ^ (e^ ® if O +1^ ® e^,) O ® if+i ® ... ® 1^ 1^. 

From [AC(e), Acq,] = 0, a = 1,2, ...,n, it follows easily that [/i, A^Cq] = 0, where 
A^e^ = (1 O 1 O A)(l (g) A)(A)e«, as defined in (|l]). Obviously [iF(/i)i,i+i, (e«)j] = 0, 
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Vj 7^ i, i + 1. Therefore we have, for all a = 1, 2, n: 



L-l 



L-1 



a) J 

i=l j=l 
L-l r i-l L 

Ca)fc ~l" + (eQ,)j-|-i 

2=1 j=l k=i+2 

L-l L-l 

= [IF{h)i^i+i, (e„)i + {eo)i+i] = E [^Wi,i+i^ {A'^ea)u+i 



(6) 



i=l 



i=l 



Let \^ be a complex vector space and R be the solution of quantum Yang-Baxter 
equation (QYBE) [g, [TT[] without spectral parameters, see e.g. |T^. Then R takes values 
in End<siV O V). The QYBE is 

R12R23R12 = R23R12R23, (7) 

where R12 = R®id, R23 = id® R and id is the identity operator on V . 

In the following we say that a ladder model having a (quantum mechanical) Hamilto- 
nian of the form 

H=Ymu+i (8) 

2 = 1 

is integrable in the sense that the operator Ti satisfies the QYBE relation (|^, i.e., 

(^)l2(^)23(^)l2 = (H)23(^)l2(^)23, (9) 

where {H)i2 = 71® id and (7^)23 = id®T-L. H is a solution of the Yang-Baxter equation 
without spectral parameters. Correspondingly the i-th complex vector space Vi now 
stands for ® Hf. After Baxterization the Hamiltonian system satisfying relation 
(^) can in principle be exactly solved by the algebraic Bethe Ansatz method, see e.g. [|I|. 

We consider ladder models with SU{2) symmetry. Let S'j, i = 1,2,3, and C be the 
generators of the algebra SU{2) and Casimir operator respectively. The coproduct of 
the algebra is given by ASi = 1 ® Si + Si ®) 1, i = 1, 2, 3. Taking into account that 
A{lF{e) = F{Aie), i = 1,2,3, j = 1,2, Ve G SU{2), the generic h is of the form 
iF(Ci,C2,C3), where 



Ci = ® 1 ® 1 ® Si + 1 ® Si ® 1 ® Si + 1 ® 1 ® Si ® Si), 

i=l 
3 

C2 = Y(^i ®1®1® Si + Si® S^®1®1 + Si® 1® Si® 1), 

i=l 
3 

C3 = ^(^2 ®1®1® Si + 1® S, ®1® Si + S, ®1® Si ®1 + 1® Si® Si ®1] 



i=l 
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In the spin-| representation of the algebra SU{2), the solutions of the QYBE (^) are 
16 X 16 matrices. For instance, it is easy to check that 

'^O ~ 108 288 *^112 H 270 *^113 H 216 *-'121 

29/^ I 90d-31/ ^ I 2d~f ^ I -54d+26/ ^ 
""108^122 H 36 <-^123 H 2~^131 H io8 ^132 UUJ 

I -108d+43/^ I -216d+80//-y i H//^ i 216d-119/^ 
+ 540 *^133 + 864 '-'211 + To8<^212 + JoS *^213 

satisfies ® for all d, f & M, where Cijk = Ci ■ Cj ■ Ck, k = 1,2, 3. 

The corresponding solution related to the 5'f/(2)-symmetric integrable ladder model 



m 



T0|] can be expressed as 

-1/ 5 11/^ 61/^ 41//^ \ I 21 

rt — — 48<-^lll — 32^112 — 30^113 ~ 48l<-^121 — ^122) + 16^123 

(11) 



,3^ 17^ , 173^ , 55^ 5^ , 131^ 

+ 4^^131 — T2*-"132 + 240^133 + 96*-"211 — 3<-^212 + ^tv213- 

Through baxterization, Tii^x) = (x — 1)7Y + 16Ii6xi6 satisfies the QYBE with spectral 
parameters: 'Hi2{x)U2z{xy)'Hi2{y) = 'H2z{y)'Hi2{xy)U2z{x) , where Hi2{-) = H{-) ® I4x4, 
'^23(') = I4x4 ® '^(■)) Inxn dcuotes the u X u identity matrix. The model can be exactly 
solved using a algebraic Bethe Ansatz method. It describes a periodic spin ladder system 
with both isotropic exchange interactions and biquadratic interactions: 

1 11 1 1 1 

H = 2 5^(2 ~^ '2^1,1 ■ Si,i+i)(- + 282,4 ■ S2,i+i) ~ 2 ^^^2 ~^ ' S2,i+i)(- + 2S2,i ■ Si,j+i' 

+7 (o + 2Si,i ■ S2,i)(- + 2Si j+i ■ 82,1+1), 
6 i=i 2 2 

where 80,4 = (a^j, cr^ j, j)/2, a^,ay,a^ are Pauli matrices. 81,, (resp. 82,4) is the spin 
operator on the first (resp. second) leg of the z-th rung of the ladder. 

It is further shown that for a more general form of (pT]), 

nil _ -45+23 a-4 6-28 , -99-3a-36-c^ , -1098-91 a-118 6-16 c 

^ - 432 "^111 + 288 '-'112 + 540 <-113 

I -369-97 g- 70 fe+50c r^ , 396+4 a+1,lb+2bc n 1 189+29 a+20fe-4c /-( 
432 '-"121 432 '-"122 H 144 ^-"123 

I 3/-y I -306-2 a-29 6-14 c /-( , 1557-71 a+172 b+l2ic r^ 
+ 4<-^131 H 216 ^-"132 H 2160 ^-"133 

I 495-a+53fe+47c /-( , -720-22 a-49 fe-43 c , 1179+91 a+118 b+16 c r^ 
+ 864 '-"211 + 432 '^212 + 433 <-213 

L-1 

with a,h,c G IR, the corresponding ladder model H' = '^'H[ j^^^ can also be exactly 



(12) 



solved by a coordinate Bethe ansatz [p!0[] . 



4 = 1 
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We consider now stochastic processes |TB[ on a ladder. Let (fi, P) be a probability 



space, with Q the finite sample space and P the probability measure defined on the a- 
algebra of all subsets of Q. For a discrete time stationary Markov chain {Xj}, i G IV, 
with underlying probability space {Q, P) and a finite state space S = {1, 2, 3, m}, there 
are transition probabilities {pajs}, a,l3 = 1,2, ...,m. The stochastic transition matrix 
P = ij>ap) has the following properties: 



m 



Pa/3 > 0, ^j9q/3 = 1, a,/? = 1,2, ...,m. (13) 

For a stationary continuous-time real- valued stochastic process, {Xt}t&]Rj^ (on the 
probability space (fi,P)), the transition semigroup P{t) = P{Xt=j\Xo = i) is generated 
by an intensity matrix Q = {qap) with the properties: 

> 0, a 7^/3, qaa = - ^qo^p, a, /3 = 1, 2, m. (14) 

The transition matrix P (resp. intensity matrix Q) defines the stochastic processes 
on a ladder. In the following we call a ladder associated with above stochastic processes, 
for instance, particles jumping randomly on the ladder, characterized by the matrices P 
and Q a Markov ladder, though geometrically it is equivalent to a chain with particular 
non-nearest neighbor interactions. If the eigenvalues and eigenstates of P resp. Q are 
known, then exact results concerning the stochastic processes, such as time- dependent 
averages and correlations, can be obtained. We say that a Markov ladder is integrable 
(resp. SU (2)-symmetric) if the eigenvalues and eigenstates of the related transition matrix 
P or intensity matrix Q can be exactly solved (resp. is SU{2) invariant). 

To every site on the i-th rung and ^-th leg of the ladder we associate states described 
by the variable r/ taking values and 1 (conventionally a vacancy at the site is associated 
with the state and an occupied state is associated with the state 1). The state space of 
this algebraic ladder is then finite and has a total of m = 2^^ states. 

[Theorem 2]. The following matrix 

1 L-i 
= 4(L-l)(18 + 4a + 46 + c) g ^^^^ 

defines a stationary discrete-time S't/(2)-symmetric integrable Markov ladder for a, b, 
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c > 0, a + 26 — 16 > 0. The operator Ti" is given by 



n" 



/ ^ 
' ai 


02 


02 


02 


03 
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04 


04 
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04 


04 


03 


04 


04 


04 


a2 


05 
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04 




Og 


05 


Og 


Og 


O4 


09 


04 


07 


Og 


03 


Og 


Og 


04 


09 


04 


02 


oe 


oe 


05 
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10+O+26, 03 




6+0+26, 


04 


= 2+0, 


05 = 


54 



(16) 



Oe = —16 + o + 26, 07 = 14 + o, Og = 8 + o, 09 = o + 2b. 7i"j_|_i acts on the i and i + 1 
rungs as defined in (^). 

L-l 

[Proof]. For the integrable ladder model with Hamiltonian if' = ^ 7^- the system 

i=l 

remains integrable if one adds to H' a constant term and multiplies H' by a constant 
factor. Moreover the spectrum of H' will not be changed if one changes the local basis of 
the rungs, i.e., the following Hamiltonian H", defined by 



H" = BH'B-\ B = ®t^Bi, 
where Bi are 4x4 non singular matrices, has the same eigenvalues as H'. 



(17) 
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It is straightforward to prove that Ti" = BTi'B , where 



B 



( 


1 





\ 


1 


1/2 


-1/2 


1 





-1/2 


-3/2 





V 


1 





-1 / 



L-1 



Therefore the Hamiltonian systems H' and H" = ^ "^^fj+i satisfy the relation (p!7| ) with 

i=l 

Bi = B, i = 1,2, ...,L. Hence H" is also S't/(2)-symmetric and integrable with the same 
spectrum as H'. 

For a + 26 > 0, as the entries of Ti." are positive, H'^^ > 0, a, f3 = 1, 2, 2^-^. From 
(|16D we also have Ea=iK/3 = 4(18 + 4a + 46 + c), V/3 = 1,2, 16. By the definition (|3|) 



Psu{2) is the transition matrix of a stationary discrete-time SU (2)-symmetric integrable 
Markov ladder. ■ 

The state space of this Markov processes associated with the stochastic matrix Psu{2) 
is S* = (1, 2, 2^^). Generally there is no closed subset C of the state space S such that 
{Psu{2))ij = for alH G C and j ^ C. In a certain parameter region of the a, 6, c there 
are nonempty closed sets other than 5* itself and the Markov ladder becomes reducible. 
However there exists no absorbing state in this Markov ladder. 

By using results in the proof of theorem 2, we have also the following integrable 
stationary continuous-time Markov ladder: 

[Theorem 3]. The matrix 

L-1 

Qsu{2)=H" -A{L-l){lS + Aa + Ah + c) = ^(H" - 4(18 + 4a + 46 + c)),,i+i (18) 

1=1 

is the intensity matrix of a stationary continuous-time Markov ladder. 

We have discussed the construction of integrable ladder models with Lie algebra sym- 
metry. It is shown that the stochastic processes correspond to the SU{2) symmetric in- 
tegrable ladder models define exactly solvable stationary discrete-time (resp. continuous- 
time) Markov ladder with transition matrices (resp. intensity matrices) which coincide 
with those of the corresponding integrable models. 

Integrable ladder models with quantum algebraic symmetry and the related Markov 
processes can be investigated in a similar way. Let e = {cq,, fa, ha}, a = 1, 2, n, be the 
Chevalley basis of a Lie algebra A with rank n. Let e' = {e^, h'^}, a = 1, 2, n, be 
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the corresponding elements of the quantum (q-deformed) Lie algebra Aq. We denote by 
Va the simple roots of the Lie algebra A. The Cartan matrix (oq/?) is then 

= ^(^a ■ r/3), da = ^{ra-ra). (19) 

The coproduct operator A' of the quantum algebra Ag is given by 

A'h'^ = h'^^l + l^h'^, (20) 

A'e^ = e'^ ® q-'^-^'- + q'^-^'- e'^, (21) 

^'fa = + (22) 

g G C, q'^°' 7^ ±1,0. In the following we use the notations A'™ and A^* defined similarly 
as in (|I|) and (1). 

[Theorem 4]. The ladder model defined by the following Hamiltonian acting in Hi ® 
Hf ® H2 (g) H2 eg) ... ® Hi H]^ is invariant under the quantum algebra Ag-. 

L-l 

Hg=J2 JF{K)i,+u (23) 

where hq = I]j=i (^ij^'i^^'/ACq^e'), Cq{e') is the Casimir operator of Aq. 
[Proof]. The generators of Aq on the ladder are given by 

L 

H'^ = A'2^-l/l^ = ^lJ®l2®...(/i;®l2®+ll®/i^)®... 1^0 12 

1=1 

L 



1=1 

L 

K = E?''"''"®-®(/a®li /a)®- 



i=l 



From [A'F{Cq{e')),A'a] = 0, Va G Aq and q^<^'-^'<^ = q^^^K q±d^h'^^ j^g^^g 
A'2/,;] = A'^e'J = K, A'V;] = 0. Therefore 



i=l 





) (e'j2 ® 


+(g'^°''-)- g 


) {q'^"'''")f 


+ (g'^"^-)| § 




E [^(^), 

i=l 


A'^(eL); 



= 0. 
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[Hq^F'^ = is obtained similarly. [Hq,H'^ = can be proved like (^. Hence Hq 
commutes with the generators of Ag for the ladder. ■ 

The Hamiltonian system (p3D is expressed by the quantum algebraic generators e' = 
{h'^,e'^, f^). Assume now that e e'(e) is an algebraic map from A to Aq (we remark 
that for rank one algebras, both classical and quantum algebraic maps can be discussed 
in terms of the two dimensional manifolds related to the algebras, see We then have 

L-l 

Hq = Y^ F{hq{e'{e)),,+,. (25) 

i=l 

In this way we obtain ladder models having quantum algebraic symmetry but expressed 
in terms of the usual Lie algebraic generators {cq} with manifest physical meaning. 
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